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Using bosonization-fermionization transformation we map the Tomonaga-Luttinger model of spin- 
less fermions with non-linear dispersion on the model of fermionic quasiparticles whose interaction 
is irrelevant in the renormalization group sense. Such mapping allows us to set up an expansion 
for the density-density propagator of the original Tomonaga-Luttinger Hamiltonian in orders of the 
(irrelevant) quasiparticle interaction. The lowest order term in such an expansion is proportional to 
the propagator for free fermions. The next term is also evaluated. The propagator found is used for 
calculation of the Coulomb drug resistivity r in a system of two capacitively coupled one-dimensional 
conductors. It is shown that r is proportional to T 2 for both free and interacting fermions. The 
marginal repulsive in-chain interaction acts to reduce r as compared to the non-interacting result. 
The correction to r due to the quasiparticle interaction is found as well. It scales as T 4 at low 
temperature. 



I. INTRODUCTION 

The bosonization has been an indispensable tool for 
one-dimensional (ID) interacting fermion studies. The 
advantage of the bosonization is that it allows to treat 
exactly the fermion interaction operator. This operator 
is marginal in the renormalization group (RG) sense, and 
therefore cannot be dealt with the help of perturbative 
approximations. 

However the bosonization becomes inconvenient when 
one needs to go beyond the marginal operators. For ex- 
ample, when calculating the Coulomb drag resistivity in 
the system of two ID wires [1, 2, 3, 4], it is necessary to 
account for the dispersion curvature v' F of the fermions 
[5] which is an irrelevant operator. In the language of 
the bosonization such operator introduces interaction be- 
tween the bosons with the coupling constant v' F . This 
destroys the solubility of the bosonized Hamiltonian. 

Moreover the Coulomb drag resistivity r is non- 
analytical in v F : 

rcc\v' F \. (1) 

This indicates that simple perturbation theory in orders 
of v' F is not applicable. 

Different methods was used to address the issue of the 
non-linear dispersion [5, 6, 7, 8, 9, 10, 11, 12, 13, 14]. 
Unfortunately, these papers either rely on numerical cal- 
culation or exact solubility or employ uncontrollable ap- 
proximations or devise methods suitable for a particular 
task at hand. No universal approach emerged from those 
works. 

A detour around the bosonization was proposed in 
[15, 16] where it was shown that a generic Tomonaga- 
Luttinger (TL) model of ID interacting spinless fermions 
may be mapped on a system of free fermionic quasipar- 
ticles with weak irrelevant (in RG sense) interactions. 



The latter approach was particularly convenient for 
evaluation of the density-density correlation function. It 
was demonstrated that this correlation function is pro- 
portional to the density-density correlation function of 
the free fermions plus small corrections due to the inter- 
actions between the quasiparticles. In Ref. [15, 16] the 
density-density retarded propagator and the den- 
sity spectral function = — 2Im£>fc w were determined 
to zeroth order in the quasiparticle interaction. 



Since the Coulomb drag resistivity is a functional of 
Bku, h is natural to apply the method of [15, 16] to 
the problem of calculating the Coulomb drag resistivity. 
This is the purpose of this paper. More specifically, we 
show below that at small temperature r aT 2 , where 
the coefficient a is a decreasing function of the in-chain 
repulsion, and that correction due to the quasiparticle 
interactions Sr vanishes quicker than T 2 : Sr oc T 4 . These 
are two main results derived below. 



This paper corrects [17] and the previous version of 
this preprint [18]. In these two publications the zero- 
temperature form of the spectral function is used for cal- 
culation of the finite-temperature drag. This mistake is 
not crucial for it does not alter neither T 2 behavior of r, 
nor T 4 behavior of Sr. It does change the 0(1) numerical 
coefficient in front of formula for r, as well as details of 
r and Sr derivations. 



The presentation is structured as follows. First in Sect. 
II we establish the mapping of the TL model on the 
quasiparticle model. Next in Sect. Ill we obtain the 
density-density propagator and the Coulomb drag resis- 
tivity to zeroth and first order in the quasiparticle inter- 
action. The conclusions of the paper are given in Sect. 
IV. Involved calculations are performed in Appendices. 
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II. THE MAPPING 

In this paper we study the model of ID fermions whose 
Hamiltonian is: 



fl 



(2) 
(3) 



Akin + H n \ + Hi nt , 

Akin = i, F /^t VVv -:4v^:), 
fl nl = v' F J rfx(:(V^)( V ^ L ):+:(V<)(V^ R ):),(4) 
= .9 J dxp R p L , (5) 



A 



where V'R.l are chiral fermionic fields corresponding to 
the right-moving (subscript 'R') and left-moving (sub- 
script 'L') fermions, pr,l =: ip R lV'r l : are cmra l fermion 
densities, and the colons denote normal ordering. The 
cutoff A is assumed for this quantum field theory. 

A non-pcrturbative approach to handle Hamiltonian 
H was proposed in Ref. [15]. There a unitary operator 
U was constructed which transforms H into the quasi- 
particlc Hamiltonian: 



H qp = fl ki „ + H nl + H[ nt + 6p{N R + N L ), 



(6) 



where A^l are the total number of right-moving (left- 
moving) fermions, flkin and H n \ have the same form as 
flkin and H n \ but with % and v' F instead of vf and v' F : 

fl kin = w F J £^:$V^ L :-:^V^ R :), (7) 
dx (: (V#)(Vi£ L ): +: (V^)(V^ R ): ) ,(8) 



H n 



v F = v F \ 1 - 



v 2ttv f / 



(9) 
(10) 



The usu al TL parameter JC is used in the last formula: 
/C = \J (2iyvf — g)/(2iiVF + g). Operators with tildes 
correspond to the quasiparticles: ip p is the quasiparticle 
field, below we use p p , which is the quasiparticle density. 

The quasiparticle interaction H[ nt in Eq.(6) is given by 
expression: 



dxp_ 



:$(Vtf p ):-:(V$)^,:](ll) 
~g> = ^(K?l*-K-^){12) 



In Eq.(ll) the summation runs over the chirality index 
p = R, L whose numerical values are p = +1 for 'R' 
and p = — 1 for 'L'. Observe that the operator H' int is 
irrelevant: its scaling dimension is equal to 3 which is 
greater than 2 [19]. (To evaluate the scaling dimension 
one has to add together the scaling dimension of two 



fermion operators, the gradient, and the electron density 
operator: 1/2 + 1/2 + 1 + 1 = 3). 

The chemical potential shift Sp in Eq.(6) is also in- 
duced by the transformation U. However since such 
shift causes nothing but additional renormalization of the 
quasiparticle dispersion parameters vp and v' F , we do not 
keep track of it below. 

The readers who are interested to learn how operator 
U is constructed should consult Ref. [15, 16]. Here we 
use bosonization-fermionization sequence to establish the 
desired equivalence between H and fl qp . That way, we 
derive the result we need with no reference to the new 
technique unfamiliar to the majority of the researchers 
in the field. 

The bosonization prescription [19, 20] expresses the ID 
chiral fermion field as an exponential of the Bose field <f> 
and its conjugate 6: 



i/>l(x) = (27ra) 



-l/2^ e iVi : [e(x)+p*(x)] ) 



(13) 



where a cx 1/A and r\ p is the Klein factor. Consequences 
of this formula are: 



Pp 



1 



2V^ 



(V$+pV9) : 



-ip: ^(VV P ): +H.c. = l -: (V$ +pV9) 2 : , 
: (W P )(VVg: ~f7 2 Pp = ^: (V$ +pV9) 3 : . 



(14) 
(15) 
(16) 



Using the above formulas we can write the bosonic form 
of H: 

H [$, 9] = H kin [*, 9] + H int [*, 9] + H nl [$, 9] ,(17) 



H k i n + Hint 



_ % f 

' 2 J 



dx 



/^(Ve^+Kr 1 :^) 2 : ,(18) 



H 



nl 



da;^:(V$+pV9) 3 :(19) 

p 



Eq.(19) was derived in Ref. [21]. 

Once the bosonic form is explicitly written we are 
ready for the second step of the derivation - rescaling 
of the bosonic fields: $ = /C" 1 / 2 *, 9 = IC 1 ' 2 ®. The 
different pieces of Eq. (17) can be expressed in terms of 
this new boson as such: 



flkin + flint - y J dx : (ve) 2 : +: (v$) * : 



(20) 



H 



nl 
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v' F J dxJ2- (/C 1/2 V<l+p/C- 1/2 V9) 3 :.(21) 



Introducing a new fermion (the quasiparticle) with the 
help of the formula: 



ftjx) = (27ra)- 1 /2 7?pe i^[e(x)+p* (a; )] ; 



(22) 



one can refermionize H. Namely, inverting Eq.(15), we 
obtain for the sum if kin + flint: 



Akin + flint — A 



kin j 



(23) 
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where -ffkin is given by Eq.(7). On the right-hand side 
of this expression the interaction term J dxp R p^ is ab- 
sent. Thus, the marginal interaction, the most trouble- 
some part of the Hamiltonian, is removed. 

The price we have to pay for the absence of the 
marginal interaction is that H n \ expressed in <l, 0, 
Eq.(21), cannot be easily fermionized. It is convenient 
to rewrite the latter equation: 

H nl = ^-i' F /dx^:(vl>+pVe) 3 : (24) 
v 

v p 

The advantage of this form is that it depends on (V$ ± 
VO) combinations only. Thus Eqs.(14), (15) and (16) 
may be immediately applied and H n \ may be fermionized: 



H n \ = H n \ + H[ nt , 



(25) 



where H n \ and H[ nt are given by Eqs. (8) and (11). 

This almost concludes the derivation of Eq.(6). What 
we lack is the dp term of H qp . To obtain this term we 
must (i) handle normal ordered expressions more accu- 
rately and (ii) take special care about the zero modes 
•Nr.l- Technically this is similar to the treatment of Ref. 
[15, 16]. Since we are not interested in the chemical po- 
tential shift, we do not address this issue here. 



III. DENSITY-DENSITY PROPAGATOR AND 
THE COULOMB DRAG RESISTIVITY 

Once the mapping of H on H qp is established, we can 
use it to calculate the density-density propagator D kuJ . 
First let us find the following density-density correlation 
function: 

K= ({ Pr (x,t) + Pl (x,t)][ Pr (0,0) + PL (0,0)}}. (26) 
In the bosonic form it equals to: 



1 



TZ= -(V$(x)e" T ' H V$(0)>. 
After the field rescaling it becomes: 



K = -(V$(x)e- rH ^V$(0)). 

7T 



(27) 



(28) 



Under fermionization 1Z transforms into: 

U = JC([p R (x, t) + p L (x, T)][p K (0, 0) + p L (0, 0)]) qp , (29) 

where the subscript 'qp' reminds that the averaging is to 
be performed with respect to the quasiparticlc Hamilto- 
nian H qp . 

One can prove through the same line of reasoning that 
D ku! =fCD ku! , (30) 



where D ku! is the retarded density-density propagator for 
the quasiparticle Hamiltonian H qp . Consequently, the 
task of finding the density-density propagator for the 
physical fermions is reduced to the task of finding the 
quasiparticle density-density propagator. The latter is 
much easier, for the quasiparticle interaction H[ nt is ir- 
relevant in the RG sense. 

As a starting point we calculate D klJ to zeroth order 



in H[ nt : 



D kuj = K [bl ku + Dl ku! ) + 0(H' int ), 
B ku = JC (j3 Rkw + B^ k Sj + 0(H[ nt ), 



(31) 
(32) 



where the chiral quasiparticlc density-density propaga- 
tors L an( i corresponding spectral densities B 1 ^ L are: 



U pkuj 



l pq-k/2 ~ n pq+k/2 



' + £pq-k/2 — Spq+k/2 + *0 ' 

B pkui = d 1 ( n pq-k/2 - n pq+k / 2 ) x 

J — A 

5(U + £p 9 -fe/2 - £pq+k/2)- 



(33) 
(34) 



Symbol f . . . stands for (2ir) 1 J_ A dq . . ., quantity e pq 
is the quasiparticlc dispersion, and n pq is the Fermi oc- 
cupation number: 



£ pq = pi>Fq + v F q 
n pq = [1 + ex P(^pq/T)]~ 1 



(35) 
(36) 



At T = the integrals in Eq. (33) and Eq. (34) can be 
easily evaluated. One finds the lowest order expression 
for D and B [15, 16]: 



D ku _ 



JC 



AlTVpk 



Bku — 



■In 



{v F k~v' F k 2 ) 2 - {uj + iQ) 2 



(v F k + v' F k 2 ) 2 - (lo + iO) 2 



(37) 



K 
2v' F k 



[tf(w 2 - (v F k-v' F k 2 ) 2 ) (38) 



-0 (lu 2 - (v F k + v' F k 2 ) 2 )] sgn cj + 0(H! nt ). 

For non-zero temperature Aristov shows [13] that 
Eq. (33) reduces to an expression with a special function. 
Eq. (34) may be calculated even for non-zero temperature 
in terms of elementary functions: 



n pku 



sinh(^) 



4|£pfc| cosh(wp ^ 



k 2 s , Sujp+v' k 2 > 



,(39) 



Suip = cj — pv F k.(A0) 



Appendix A provides the details of the calculation. 
Eq. (39) can be used to find B kuJ , Eq. (32). 

Once the spectral function is found the Coulomb drag 
resistivity may be evaluated with its help. Before pro- 
ceeding with such calculations, let us briefly explain what 
Coulomb drag is. 
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In the Coulomb drag experiment two parallel ID wires 
(subscript i = 1,2) of length L with Hamiltonians H, 
Eq.(2), are coupled capacitively with the Hamiltonian 
He = 9c J dxpip2- Because of this coupling, electrical 
current / in one of the wires induces potential drop V 
across the other wire. The proportionality coefficient 
between V and / is called the Coulomb drag resistiv- 
ity r = V/IL. It characterizes the pulling force, which 
the fermions in the current-carrying wire exert on the 
fcrmions in the other wire. 

The experimental physics of the Coulomb drag is quite 
rich: the observed values of r could be cither positive or 
negative, and show dependence on temperature, spacial 
inhomogencity, and applied magnetic field. 

In general, two mechanisms are discussed in the the- 
oretical literature [3]. According to one mechanism, the 
drag occurs because Wigner crystal-like correlations in 
both wires lock against each other. In the RG lan- 
guage this corresponds to the relevance of the inter- wire 
backscattering interaction. Such mechanism works best 
at zero temperature but quickly deteriorates at T > 0. 
We do not study this mechanism in our paper. 

The second mechanism is insensitive to the inter-wire 
backscattering. Instead it relies upon interaction of the 
smooth components of the electron densities. It is more 
resilient towards temperature but the resultant value of 
r is proportional to \v' F \. Because of this, it is impos- 
sible to study the second mechanism within the well- 
established framework of the one-dimensional bosoniza- 
tion. The purpose of this paper is to provide a reliable 
approach overcoming this difficulty. 

We start with the following formula for r (Eq.(7) of 
Rcf. [5]) which is valid when the inter- wire backscatter- 
ing can be neglected: 



9l 



167r 3 n 2 T 



dk 



+oc \ k 2 B 1 ^B. 
aw 



2.ku> 



sh 2 (w/2T) 



(41) 



Here n is the electron density. The spectral density 
A(k,oj) of Rcf. [5] relates to our spectral density B^ 
as: Bkij = —2A(k,u). In that reference the notation U12 
is used for the inter- wire coupling constant gc- Below 
we study the case of identical wires: B\ = B 2 = B. We 
assume that \vL\ is small: 



|u F |A < v F . 



(42) 



Since we know Bku> to zeroth order in g' 
forward to find r with the same accuracy, 
our calculation with the following observation. The func- 
tion: 



it is straight- 
Let us begin 



bpku. 



a pkuj 



shm(#) 



(43) 



\oj -pv F k\ < v' F kT/v F < v' F T 2 /v F . 



(45) 



Outside of this region b p ku, is exponentially small. 

Since b-^ku is exponentially small when k and oj are pos- 
itive, for the purpose of calculating the integral Eq. (41) 
we may further approximate: B w JCB^. Substituting 
this in Eq. (41) we obtain: 



-AT), (46) 



256tt 3 (v f ) 2 ti 2 T 
r+00 r+00 
I(T) = 16(5 F ) 2 / dk / duk 2 bl kuJ (47) 
Jo Jo 



dkduj 



cosh 2 (v F - 



-) cosh 2 (%>- 



iv'kT 



The integral I(T) is evaluated in Appendix B. It is shown 
that, if temperature is low: T < 6j?, where renormalized 
Fermi energy is: 



then 



This gives us: 



32|5 F |T 3 



„~,3 



+ o(T 3 ). 



r = 



87T 3 W' 



j ,3 L ,/ 



,-T 2 + o(T 2 ). 



(48) 



(49) 



(50) 



F I F 



The electron zero-temperature density, which enters Eq. 
(50), is not an independent quantity. It can be expressed 
in terms of the dispersion parameters: n = v F / (2ir\v F \) . 

Furthermore, it is possible to show with the help of 
Eqs. (9) and (10) that v F = v F + 0(g 2 ) and v' F = v' F + 
0(g 2 ). Since our accuracy does not allow us to keep 
0{g 2 ) terms, we can replace v F by v F and v' F by v' F in 
Eq.(50). Thus, it is true: 



aT 2 , 
9 2 C IC 2 \ 



2ttv f 



+ 0(g 2 ). 



(51) 
(52) 



The interaction enters the expression for a through JC 2 = 
1 — g/(nv F ) + 0(g 2 ). We see that the repulsive in-chain 
interaction (g > 0) acts to reduce r. 

Eq.(51) may be cast in the following form: 



ci/C 2 
lo 



(53) 



2\v' F k\ [cosh(|g^)+cosh(f^) 

localized mostly at small momenta and near the 'light 
cone': 



|fc| <T/v F , 



(44) 



where l 1 and e F are defined in Ref. [5]. They are: 

e F + 0(g 2 ), (54) 

(55) 



1 

To 



9c 
2irv F 



£f 
2 



9l 



8^ 3 \v' F \v F 
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The numerical coefficient c\ is equal to 7r 2 /4. For free 
fermions (JC = 1) the formula identical to Eq. (53), with 
the same value for c\ is derived in [5] . 

Furthermore, that paper establishes that r oc T 2 for 
both non-interacting fermions and for exactly soluble 
Calogero-Sutherland model. Our Eq.(50) proves that T 2 
dependence holds for a generic interacting model as well. 

In Rcfs. [17, 18], where, mistakenly, the zero- 
temperature form of B is used in the integral for r, differ- 
ent value of the coefficient C\ is found. Present derivation 
corrects that error. 

Eq.(51) and Eq. (52) account for effects of marginal 
interaction g. In addition to that, our method allows us 
to evaluate the correction Sr to r due to the quasiparticle 
interaction g'. 

To find Sr we calculate the lowest order correction to 
D koj . Such correction appears in the first order in g': 
since H[ nt couples the right-moving and the left-moving 
quasiparticles, the expectation value (pr(x, t)pl(0, 0)) 
taken with respect to -ff qp , Eq.(6), is no longer zero 
but instead 0(g'). Thus, full Matsubara propagator 
V = V° + SV, where: 



SV(x,t) = K5V{x,t) 



(56) 



K / dr' 



({pn(x,r)HUT')p L (0,0)))o 



+ ((p L (x,T)H! nt (r')M0,0)))o 



Here the symbol ((. . .))o stands for time-ordered averag- 
ing with respect to the non-interacting (g' = 0) quasipar- 
ticle Hamiltonian. Matsubara propagators are denoted 
by calligraphic letters (e.g., T>, V), retarded propagators 
are denoted by italic letters (e.g., D, P). 

Feynman diagram, which describes ST>, is shown on 
Fig.l. The wavy interaction line is to be identified with 
g'(k\ + k 2 — pi — P2), h, Pi are the fermion momenta. 

When evaluating 5T> we note that in the non- 
interacting quasiparticle Hamiltonian the left-moving 
and the right-moving quasiparticles are decoupled from 
each other. Consequently, the object ((/9L,R^q' n t/?L,R.))o is 
split into products of left-only and right-only expectation 
values. 

As a result of this derivation one can show that the 
lowest order correction to the Matsubara propagator due 
to the quasiparticle interaction is: 



(57) 



where t> pkuJ is the chiral non-interacting Matsubara prop- 
agator, the propagator V p is defined as: 

V p (x,r) = -ip((p p (x,T) (58) 

x {:$((), 0)[V^ p (0, 0)]: -: [V$(0, 0)^(0, 0): }»„. 

It equals to: 



- 1 , (iu-pvFk) ^ 



■2kK- ' Pv' F k ~ pkuJ ' 



(59) 



To find the correction to the spectral function SB 
it is necessary to perform the analytic continuation in 
Eq.(57): 5D ku = -2g'JCJ2 p D°_ pkul P pklv , where the re- 
tarded quantity P pkuJ = V p ku\iw^w+i0 is equal to: 



pku 



1 , (uj - pv F k + iO) ^ 



2irv 



pv' F k 



pkuj • 



We obtain for SB: 
g'K 



SB k „ = ^tJ2 B-vku Unv F ReD° pku! + 1 



(60) 



(61) 



This is the lowest order correction to the spectral density 
due to the quasiparticle interactions. 

The square of the spectral density, which enters the 
equation for the Coulomb drag, is: 



B 2 



(1CBL) 2 (62) 



2~g>lC 2 



+ 



=T~ E B° p ^B°_ pku (4™ F Re^ + l) . 



We already pointed out the fact that b\, ku] is exponen- 
tially small for positive k and u>. Therefore, in the region 
of the integration we may neglect the terms proportional 
to B°B R and (S°) 2 : 



Bl « (ICB° Rk J 2 + ^(B Rk J 2 (Anv F ReDl ku , + l)(63) 



The correction to the drag is: 



g 2 c ~g'ic 2 

8n 2 v F v 2 T 



Sr = 



dkdu) ( 4irv F ReD° 



hku 



cosh(||^) + cosh(^) 



2v' F kT 



, (64) 



Since the denominator in this integral is very small ev- 
erywhere except near u> = v F k, we need to know the be- 
havior of RcI?l in this region. It is proven in Appendix C 
that: 



4™ F ReD° fcw + 1 



Slot- 



2v F k 



1 

12 



Vpk 
v F 



(65) 



1 /^R. X 2 

4 \ v F k 



^T 2 

-F 



• 1 - -2iir +0 < r3 >- 



where it is assumed that v F k = 0(T), and vpSu^/v^k = 
0(T). 

At T = the above expansion becomes: 



47re F Re£>° fc , 



1 



(Su; p ) 2 
2i F k 4vU 2 



pSu>p 



(v F k) 2 
12v 2 



(66) 



This expansion is derived in [17, 18], see Eq. (51) of these 
references. 

Finally, placing Eq. (65) into Eq. (64) one finds 



Sr cx T 4 . 



(67) 
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This result may be obtained without performing actual 
integration. It is enough to use the dimensional analysis: 



Airv F BcD° pkuJ + 1 = 0(T 2 ), 
dk = 0(T), 
duj = 0(Suj R ) =0(T 2 ). 



(68) 
(69) 
(70) 



From these Eq. (67) follows. 

Consequently, at low temperature the interaction cor- 
rection Sr to the Coulomb drag resistivity r vanishes 
quicker than T 2 . Thus, the non-interacting quasiparticle 
result, Eq.(50), suffices to capture the leading behavior 
of r at T -> 0. 



IV. CONCLUSIONS 

Using bosonization-fcrmionization trick we mapped 
the Tomonaga-Luttinger Hamiltonian with non-linear 
dispersion on the Hamiltonian of the quasiparticlcs with 
irrelevant interaction. This mapping allows us to eval- 
uate the density-density propagator of the Tomonaga- 
Luttinger model with non-linear dispersion. The propa- 
gator itself was used to calculate the temperature depen- 
dence of the Coulomb drag resistivity r. It was estab- 
lished that r oc T 2 at low T for both interacting and free 
fermions. The irrelevant quasiparticle interaction intro- 
duces additional correction which vanishes as T 4 . 
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Substituting q = (u) — pvFk)/2v' F k, as required by the 
delta-function, into expressions for i pq +k/2 we obtain: 

, 5u v + vLk 2 
, + fc/2 = ^p^, (A3) 

(Suj p + v' F k 2 ) (Suj p + 2pv F k + v' F k 2 ) 
£ m+ k/2- ^ , (A4) 

Sujp = uj — pv-pk. (A5) 

One can check that the right-hand side of Eq. (A4) van- 
ishes both at 8ui p = —v' F k 2 and at 8oj p — 2v F k — v' F k 2 . 
Of these two roots only the first is physical, the second 
is spurious. The spurious root appears due to unjustified 
extrapolation of Eq. (35) to large momenta ~ vp/v' F . In- 
deed, using Eq. (A3), we find that the momentum q + k/2 
is close to zero for the physical root and close to — pvy / v' F 
for spurious root. Since our treatment is valid near the 
Fermi points only we replace: 



Sujp + 2pv F k + v' F k 2 w 2pv F k. 



(A6) 



In the above equation we neglected the terms, which are 
quadratic in T, see Eq. (45). Consequently: 



£p9+fc/2 ~ 2^p ( p Vp >' 
With the help of the same arguments one derives: 

£pq - k/2 2^p ( p_% >' 



(A7) 



(A8) 



The expressions Eq. (A7) and Eq. (A8) allow us to per- 
form integration over q in Eq. (A2): 



°pkuj 



sinh(#) 



4\v F k\ cosh(% ^'fcT fc2 ) cosh(% ' 5 " 4 i ^fcT fc2 ) 



(A9) 



This is the quasiparticle spectral function for T > 0. 



APPENDIX A: SOME DETAILS OF B° kw 
CALCULATION 



APPENDIX B: CALCULATION OF I{T) 



Here we fill the gap between Eq. (34) and Eq. (39). 
Observe that: 



sinh 



ipq-k/2 — n pq+ k/2 



£ pq + k/2 £ pq-k/2 



2 cosh ( £ -^a \ cosh (££^ 



<A1) 



The integral I(T) is defined as: 

f+oo r+co 



r+oo r-t-oo 

I(T) = dk deux (Bl) 

Jo Jo 



1 



COSh (V F 4c , k7 r )cOsh (V F 4 -, fc / ) 



The delta-function in Eq. (34) enforces i pq+k / 2 — 
£ P q-k/2 = v- Therefore: 

= I ^ Sinh(# ^ ( ::^I 2 ^ (A2) 



pkui 

sinh 



2 cosh(2q^) cosh(^^) 



dq 



2|«F fc l J 2cosh(%^)cosh(^Zi) 



In the above integral we introduce the new variable: 



n 



VF 



V F Jo 



W F \kT 



I= mir d k 1 



(oj-v F k- \v' F \k 2 ),(B2) 
kdfl 



n cosh 2 ft cosh 2 (ft + 3gfc ) 



ftn = UF 



v F v F k 
415' IT + 4T 



(B4) 



7 



If temperature is low: T <C ep, where 



From which one gets: 



4^' 



(B5) 



then f2o ^> 1. In such a situation the integrand is expo- 
nentially small for Q, < — Oq: 

cosh -2 ftcosh~ 2 (f2 + vpk/2T) < cosh -2 fl (B6) 
< 4cxp(-2O ) < 1. 

Therefore, we may extend the integration interval: 

r>+oc /> + oo 



/-l-oo /-l-oo 
/ dfi.... 



The integral J can be expressed as: 



7 = 



/ fcdfc / 

JO ./-< 



Jo J-00 COsll ttcosll ^) 

16|6 F |T 3 /"+ 00 dQ f +oc QdQ 



(B7) 



<B8) 



cosh 7 cosh (0 + Q) 



Integrating over Q by parts we derive: 



I(T) 



16|£ F |T 3 



„~,3 



-4^An). (B9) 

cosh Si 



/■+oo 

/(O) = / [1 - tanh(fi + Q)} dQ. (BIO) 
Jo 

Since /(O) — > when f2 — > +oo, we can integrate by 
parts over f2 as follows: 



16 £r T 3 f + 

J(T) = ^ — / (tanhft + l)f'(il)d&l, (Bll) 

W F J-oo 



Therefore: 



F 

+00 



o cosli (ft + Q) 



= tanhft - 1. (B12) 



+ 00 



J(T) = 16| !^ |T3 / (1 - tanh 2 Q)dSl. (B13) 

W F J-oo 

The last integral can be calculated easily, and 7(T) is: 

32|5 F |T 3 



I(T) 



?,3 



(B14) 



This concludes our derivation of I(T). 



APPENDIX C: EXPANSION OF ReD L 

In this Appendix we evaluate RcDl kLU when lo ss v F k. 
We start with the following expression: 



ReDl kui = Rc / ^ g - fc /2- ^ +fc /2 (C1) 



p.nO 1 p „ / n Lg-fc/2 ~ ^Lq+fc/2 . , 

In the last expression we integrate by parts to obtain: 



1 



In 



y ~ 2 2lO; 



« — - — 

y 2 25TI 



(in 



L 9 



(C3) 



Since we want to derive Re-D^fe^ near the point <5wr ~ 0, 
it is convenient to substitute Sujl with the equal quantity 
5wr + 2v F k. After such substitution and some transfor- 
mations the expression for ReD^ kbJ becomes: 



In 



VF 2{?f 2vpk 



i?F 2{jp 2vpk 



dn Lq 
dq 



(C4) 



Logarithm in this formula may be expanded in orders of 
small parameter T/e F . Indeed, it is true that = 
O(T) and i F \5cu R /i F k\ = 0(T). Observe also: 

J dq J \i F i F J cfel + cxp(e/T) V ; 

v' F f £ 2 _ tt 2 v' f T 2 



cosh 2 (e/2T) 3v F 
where we used q w — £l 9 /% + ^f^l 9 /^f an< ^ the identity: 



x dx 



cosh a; 6 
Using similar procedure one proves: 
dn La tt 2 T 2 



/ 



q dq- 



dq 



(C6) 



(C7) 



Therefore, both q and q 2 , upon integration, contribute 
terms of order T 2 . Thus: 



RcD 



1 (v' F k 



o =_J_ / 



dn Lq 
,.: dq 



l + ^-^(C8) 



•- / "^ 2 1 /<W ' - 



Up 



JR 

4 V w F fc 



0(T 3 ). 



12 V % 

Since j q dn^ q /dq = l/(2n), Eq. (C8) may be written as 



4irv F ReD° + 1 = 



1 (v' F k 



2v F k 12 V v F 



(C9) 



1 /(W R \ 2 tt 2 fT\ 2 ^. m3 , 
"lid) "24 W + ° (T ) - 



This is the expansion we need. 



8 



R k L 

k 2 p 2 



FIG. 1: Feynman diagram corresponding to the lowest order 
interaction correction to the density-density propagator. 
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